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Abstract —Wireless indoor positioning systems have been 
shown to be very useful in many applications and have 
been the subject of a considerable amount of research, 
mainly concerning the two-dimensional (2D) case. 
However, in many practical situations it is necessary to 
determine the three-dimensional (3D) coordinates of an 
object or user. In this paper, a hybrid algorithm for 
implementation in a 3D indoor positioning system is 
proposed. This algorithm is implemented by using a 
fingerprinting technique based on both the k-means and 
naive Bayes methods, and uses the received signal 
strength (RSS) as an input parameter. In addition, a 
comparison of the main algorithms discussed in previous 
research papers and the proposed algorithm is presented. 
Indoor positioning experiments were conducted in a 
typical building with two floors (180m 2 ) and four access 
points (APs). 

The proposed algorithm exhibited a better performance 
than that of other algorithms, with a mean error around 
1.80m. 

Keywords — 3D Indoor positioning, Fingerprint, Naive 
Bayes, k-means, Cost 231 multi-wal. 

I. INTRODUCTION 

Technological development has made the location of both 
people and objects in outdoor environments possible by 
means of tools such as the Global Positioning System 
(GPS). However, these tools are not efficient for 
handling indoor locations because the variability of these 
environments is high as compared that of outdoor 
environments. This is because the variability depends not 
only on the antenna’s characteristics, but also on the type 
of construction and internal structures, such as the walls, 
floors, and partitions walls. 

Because of the numerous applications of indoor 
positioning systems (IPS), including utilization in 
emergency systems, localization of mobile robots, and 
navigation assistance in malls, schools, university 
campuses, airports, and hospitals, several studies have 


been conducted on location estimation in this type of 
environment. 

A survey of the main implementation techniques shows 
that they can be classified into three typical location 
estimation schemes: triangulation, scene analysis 
(fingerprinting technique), and proximity. These schemes 
are discussed in [1] and [2], whereas in [3] and [4] the 
application of artificial neural networks (ANNs) to the 
problem was proposed. The k-nearest neighbours (k-NN) 
algorithm and the Bayes method were addressed in [5], 
[6], [7], [8] and [9]. 

In this paper, we propose a three-dimensional (3D) 
solution to this problem through a hybrid algorithm. This 
algorithm is implemented by using a fingerprint technique 
based on both the k-means and naive Bayes methods, and 
uses the received signal strength (RSS) as an input 
parameter. Furthermore, we compare the proposed 
method with the ANN multilayer perceptron (MFP) and 
radial basis function (RBF) algorithms, as well as with the 
k-NN and Bayes (histogram and kernel) methods and the 
method of triangulation (lateration), which is one of the 
most commonly used methods in the context of outdoor 
locations. 

The experiments were conducted in two distinct 
scenarios, both in a typical build with a total area of 
180m 2 (first floor 131m 2 and second floor 49m 2 ) and four 
APs. In the first scenario, actual measurements were 
performed and the RSS values obtained in pre-defined 
regions. In the second scenario, the signal behavior was 
simulated from the Cost 231 multi-wall (MWM) model. 
The rest of this paper is structured as follows: Section II 
presents an overview of the main techniques and 
technologies related to indoor positioning, proposed 
algorithms are discussed in section III, section IV presents 
the main discussions and computational results and 
finally, section V summarizes the results and presents 
suggestions for future work. 
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II. INDOOR POSITIONING SYSTEMS 

2.1 METHODS FOR INDOOR POSITIONING 
SYSTEM 

In literature, the indoor positioning system can be 
classified into three main categories: triangulation 
(method based on geometric properties), fingerprint 
(signal data collection or pattern recognition 
technique) and proximity. A comparison of these 
methods is presented in [1]. 

2.1.1 TRIANGULATION 

The Triangulation method uses geometric properties 
of triangles to estimate the target location [1] and [2]. 
This method is classified into two categories: 

• Lateration: estimates user location by calculating 
the distance between the mobile unit and a set of 
APs with known coordinates. For a 3D 
positioning, it is necessary at least four APs. 
Centered on each of the APs, a sphere with radius 
r, being the location of the user defined by the 
intersection of such spheres. Figure 1 shows this 
concept. 



3D INDOOR 
POSITIONING 



Fig. 1: Indoor Positioning System Based on Lateration 
• Angulation: The location of the target can be 
estimated by the intersection of several pairs of 
lines with respect to the direction of the angles, 
each formed by the circular radius from a set of 
APs. Figure 2 shows this concept. 



Fig. 2: Indoor Positioning System Based on Angulation 
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The most common methods used to estimate distance in 
the triangulation method are [1]: 

• Time of Arrival (TOA): The distance between the 
mobile unit and the unit of measurement is directly 
proportional to the propagation time. For 3D 
positioning, TOA measurements must be 
performed in relation to the signals of at least four 
APs. 

• Time Difference of Arrival (TDOA): It consists of 
a variation of the TOA, however, the estimation of 
the mobile unit is obtained from multiple 
measurements in order to obtain the difference 
between the arrival times of the signal to the 
receiver. 

• Received Signal Strength (RSS): This technique is 
based on attenuation of the signal between 
multiple transmitters and the receiver. The 
estimation of the difference between transmitted 
and received RSS is determined through empirical 
and theoretical propagation models. Some of these 
models are discussed in [10]. 

A comparison of the methods in question, including 
parameters such as cost and accuracy, is presented in [11]. 

2.1.2 FINGERPRINT 


Fingerprint is a technique based on pattern recognition 
and involves the division of the location system into two 
phases, offline and online phases. In the offline phase, the 
vectors of received signal strength are collected the of all 
the detected Wi-Fi signals from different access points at 
many known locations (reference points) [8]. The online 
phase consists of reading an RSS value, and by using a 
classification algorithm, this RSS value is compared with 
the values stored in the database during the offline phase, 
thus obtaining the location. The Figure 3 summarizes this 
procedure. 



OFFLINE PHASE 


ONLINE PHASE 

RSS READ A ‘ 



POSITIONING INDOOR 
2D OR 3D 


Fig. 3: Indoor Positioning System Based on Fingerprint 


2.1.3 PROXIMITY 

Among the three methods for indoor location, the 
proximity method is considered the simplest. Location 
systems based on this method do not provide a coordinate 
as a user location estimate. In this case, the indoor 
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position is given sectorally, that is, the system returns the 
room in which the user can be at a given time. 
Technologies used by this method include radio 
frequency identification and Bluetooth [1]. 

2.2 TECHNOLOGIES FOR INDOOR POSITIONING 

2.2.1 Radio Frequency Identification - RFID 

It consists of a technology that uses radio waves to 
identify people or objects. Identification can be obtained 
by reading information stored in a tag attached to an 
object. Some works that use this technology include [13], 
[14] and [15]. 

2.2.2 Wi-Fi 

It is a technology that uses radiofrequency to transmit 
data from devices based on one of the 802.11 standards 
developed by the IEEE. These devices include for 
example: Notebooks, cell phones, cameras and TV's, 
which makes this technology more suitable for 
implementing an IPS [16]. 

2.2.3 Bluetooth 

It is often classified as a wireless personal area networks 
(PANs) technology, with transmission power and range 
much smaller than a Wi-Fi network, occupying a limited 
space around the user (usually 10 meters). It is used to 
connect wireless devices within a short distance, such as: 
Cellular phones, TV's, notebooks and stereos. For more 
detailed discussions on the application of this technology 
to the problem in question, see [17], [18], [19], [20] and 
[ 21 ]. 

III. PROPOSED SOLUTION AND GENERAL 
DESCRIPTION OF THE SYSTEM 

3.1 METHODS FOR INDOOR POSITIONING 
SYSTEM 

The proposed algorithm is based on Bayesian inference 
and use the fingerprint technique as the basic structure for 
estimating the location. The algorithm consists of a 
combination of the k-means and the naive Bayes, a 
simplified version of the Bayes theorem discussed in the 
subsection 3.3. The Bayes theorem is given by [22]: 

IW.)- (1) 

E P(y\x t )Pln) 

t I 

The central idea of the proposed method is to determine 
the most probable position of the target based on the RSS 
vector measured in relation to all n APs. This procedure is 
performed through the following rule: 

Let©!, S 2 ,..., Sjv be a collection of subsets, called 
sectors or rooms, where each is a candidate for target 
location and consists of n RP’s where= 0 and 
U i =1 &i = % represents the indoor space for the set of 
sectors and s = {RSSi, RSS%-- -,RSSn} the RSS vector 
observed during the online phase. Select S*, if 


P(Sj|s)>P(S ; |s), for i,j = 1,2,..., N and i ^ j. This results 
in the following model: 


P(6,|s} 


Pis|6 t )P(<Bj) 

E P(H|6 t )P(6e) 


( 2 ) 


{ 1 

When the sector that represents maximum probability is 
selected, the 3D location of the target is defined from the 
centroid of the RP’s belonging to the sector in question. 


3.2 ESTIMATION OF THE LIKELIHOOD FUNCTION 


Let x = {x\ ,X 2 , • • • an RSS sample obtained from an 
unknown probability distribution with density fx(x), in the 
off-line phase of the location system, our objective is to 
estimate the density in question. There are several 
methods for this estimate. The most common are 
presented below: 

• Histogram: this method subdivides the space of 
fxix) with support [0,1] in n bins of (Sectors oM) 
equal sizes, given by: 


Bi=[0, £].• •• .ftf-pn 1 - 1] (3) 


To estimate P(sl6;), we obtain the frequency 
distribution of the signals in the environment. The 
goal is to check the frequency of each RSS interval 
for all locations. 


Density estimation Kernel (EDK): In probability 
and statistics, EDK, is a non-parametric method for 
estimating a probability density function fx(x) given 
by: [9]. 




( 4 ) 


Where rj is the sample size, h is a non-negative data 
smoothing parameter called bandwidth. Similar to 
the histogram method, the kernel estimation 
establishes a function for representing the 
probability distribution, with the difference that the 
data is not to be allocated in discrete regions (bins), 
instead a continuous function is defined. K(*) is the 
kernel function in which: 


/. 


K\{x)dx = 1 


( 5 ) 


A widely used function in this case is the Gaussian 
kernel [9].Thus, the estimate for the term P(slS;), 
can be modeled through a normal distribution is 
given by: 


fA*) = -k£ifa**P 






( 6 ) 


3.3 PROPOSED SOLUTION (KMEANS-BAYES) 
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The implementation of an IPS through this algorithm 
divides the off-line phase into two steps. The first 
step starts with reading the RSS vectors in n RP’s 
with predetermined 3D positions that make up the 
vector r. For each RP, measurements are taken for all 
APs. In addition, these measurements are performed 
at different heights and with the mobile device 
directed to the north, south, east and west. This 
results in the following matrix: 

/RPi(-i) RPi(z i) ••• RP.\{z i)\ 

. m\ I RP> .. •• /.*/'•, ... 

A ■ (7) 

\RP[(z.\) RP 2 (zx) ••• RPx(z.x)/ 

This procedure aims to make the database robust, thus 
maximizing the accuracy of the algorithm. When all 
measurements are completed, in the second step, the k- 
means algorithm partitions the indoor space into P 
clusters that represent the sectors in which the RP’s will 
be allocated. For each of the clusters a set of centroids 
with Cartesian coordinates randomly chosen in the first 
iteration is defined. From the second iteration, the 
centroids are obtained based on the mean coordinate of 
the RP’s. In addition, a value is defined for P is defined. 
Then, the distance between each RP and the centroids is 
calculated by assigning each RP to the cluster that has the 
nearest centroid. This procedure is repeated until there is 
no RP change between the clusters. 

This procedure is performed by minimizing an 
objective function, called the squared error function, 
which is defined as 


V A’ 

• ; =EZ H r > 



( 8 ) 


i-l >=1 

Where 11 r^ 11 is a 3D measure of the distance 

between each reference points and the centroid.In general, 
Euclidean distance is applied in (8). This, together with 
other measures of similarity, including Mahalanobis, 
Minkowski, and cosine distances, are discussed in [23]. 

When the clustering process is completed, there 
have the following configuration: 

• Each RP has a definite 3D coordinate P = (x,y,z )and 
an RSS vector related to this coordinate. 


• are assigned probabilities for each of the clusters, 
according to the knowledge of the frequency of 
users in each of these clusters or evenly, that is, 
considering that POP*)** = 1,*** ,N are equiprobable. 
After this procedure, the offline phase is complete. 

In the online phase the vector s = 
{RSS\,RSS%--,RSSn} is read and the target position is 
estimated based on equation (2), resulting in: 


P\% |s) 


P(s) 


P( RSS U RSS,. • • • . RSSs % )P('Vi : 
Pi S| 


(9) 


Using the general multiplication rule [24] in (9), 
results in: 


P I RSSi . RSS 2 , • • • . RSSs |«p.) = 

P{RSSi | RSSi. • • • • RSSs . %) x 
K P<RSS.\RS .<■.■■■ ,RSSs,% . ( 10 ) 

X Pi RSSs _, I RSSs ,Li). P( RSSs flft) 

In order to reduce the computational cost of system 
incurred by equation (9) when it is necessary to acquire a 
large volume of data, we will consider that that the 
probability of each attribute of vector s are independent of 
each other given the locations. This means that 
P(RSSi,RSS 2 ,--- ,RSSn\^i) can be simplified: 

jV 

nw,i do 

; i 


Thus, the target position given s is simplified to: 
l’ %\\s) a P\%. RSSi. - - RSS,\ ; 
x P( RSS ]. • ■ • , RSS,\ %) 
a P(%) x P[RSSi\%) x x P{RSS s \%) 

X 


a Pi%\ U/’[«!$') 

t=l 


( 12 ) 


Algorithm 1 Hybrid algorithm: k-means and naive Hayes far 
t hree-dimmsinnal indoor positioning in wireless local arrn networks 


Input; r - IRR.,RP>. . RP,\) e R n that iepret*>ut« the vector 
of rHerence points. 

Number of clusters ^3 reprreenting the number of ©ectors 
A wt»n) vector a — (RSS\, RSSi, • . RS$s), 

Output; A 3 D posit urn from the centrucd coordinate of the HP's 
t Define a set of centroids t*j, cj, , op Define a set of centroids 
with random coordinates for each cluster. 

7 reju'Ht 

x Calculate distances between the RP’s and the centroids ruing 
equation (h) ami assign each point to the cluatei with the 
client centfold. 

i: Recompute the centroid of each duster tuord oil the coordinates 

of tile FCP's* 

v until converged (there is no movement of data between the dusters) 
« for i->l to N do 

V. Compute PCPtjsJ Itosed oo equation ( 10 ) 

>: end for 

Allocate Use target to the cluster that results m the maximum 
probability. 

in Cct the lil 1 ’* cmtrokd. 

it Return 3 D portion Lug tu*Mvl on the oxirdinaU*? obtained in step in. 


Equation (12) is a simplified version of the Bayes 
theorem, known as the Naive Bayes. More details on this 
classifier can be obtained in [22]. After the cluster that 
result in maximum probability has been identified, the 
target position is estimated from the centroid coordinate 
of the RP’s. Algorithm 1 and figure 4 summarize the 
proposed method. 
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IV. Experimental Results and Discussion 
4.1 INTRODUCTION 

The experiments were conducted in two distinct 
scenarios, both in a typical build with a total area of 
180m2 (first floor 131m2 and second floor 49m2) and 
four APs. In the first scenario, actual measurements were 
performed and the RSS values obtained in pre-defined 


regions. Figure 5 shows the floor plan with for the two 
floors used in experiments. In order to make the 
developed IPS more realistic, all experiments were 
performed in the environment with people performing 
their normal activities. 
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3500,CO ►35CO.QO ► 2059,10 ► 



9321.04 ► 

Fig. 5: Floor Plan used in experiments and simulations 


In the second scenario, the same build was considered, 
taking into account all the characteristics of this 
environment (for example, the quantity, position and 
transmit power of the Aps), with the difference that in the 
scenario in question no real measurements were Instead, 
the Cost 231 multi-wall model (MWM) was used to 
simulate the behavior of the signal in this residence. This 
model considers the signal attenuation as the loss in free 
space added to the loss resulting from the numbers of 
walls and floors penetrated between the transmitter (T x ) 
and receiver (R x ). The model in question is given by [10]. 

' & 

£ - Llf + L f + 2_, button + kj Lf 

(13) 


where 


• LFS = free space loss between transmitter and 
receiver; 

• Lc = constant loss; 

• kwj = number of penetrated walls of type; 

• kf= number of penetrated floors, 

• Lwi = loss of wall type; 

• Lf - loss between adjacent floors, 

• h- empirical parameter; 

• l- number of wall types. 


The total loss related to the floors is a non-linear function 
of the number of penetratedfloors[10]. This feature is 
taken into account by the introduction of an empirical 
factor h. For practical reasons, the number of distinct 
walls should be kept low because the difference between 


the various wall types in relation to the attenuation is 
small and its meaning in the model is not clear [10]. The 
types of walls in question are displayed in the table 1 
[10], while the typical values of Lwl, Lw2, Lf, and b are 
shownin the table 2 [10]. 


Table. 1: Types of walls used in the Cost 231 MWM 


Wall lypc 

Description 

l .iytn wull </^|) 

A wnll that it not *>eannp toad: c.g. plasterboard. panicle board ir 
ihin MO cm). Iijjtit concrete wall 

Heavy wall (/. w2 l 

A lujcJ-neJnng wall or ocber th>ci I 10 cm| wall. 

•tniik* i>( c £ concrete or brick 


Table. 2: Values for Lwl Lw 2 , Lw2, Lfeb in dB 


Environment 

Multi-wall model (MWM) 

Lwl [tlBl 

uzijbi 

L f |«JB| 

b 

Dense 

one floor 
two floors 
multi floor 

3.4 

69 

183 

046 

Open 

14 

69 

183 

0.4* 9) 

Large 

3.4 

6.9 

183 

0.46 

Corridor 

3.4 

69 

183 

0 46 


4.2 PARAMETERS OF ALGORITHMS 
As algorithms for comparison, the main solutions 
addressed in the literature were chosen. The following 
parameters were used: 

• Triangulation: Let D = (x, y, z) be the coordinate 
representing the location of the target, given by 
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the intersection of four spheres centered on four this case, the IPS was developed through equation (2). 


APs with known coordinates given by (xi, yi,zi), 
(X 2 , y 2 ,Z 2 ), (X 3 , y 3 ,Z 3 ) and (x 4 , y 4 ,z 4 )whose 
distances to the target are given by di 2 with i = 
1,2,..., 4, estimated by the attenuation of RSS 
,we obtain the location of the target from the 
following relation: 

(x - Xi) 2 + (1/ - j/i)- 4- (z - z\ }* = dj 

j t 1 

,2 i < V’ l2 (14) 

[x a 

(x - X 4) 2 + (j/ - y. i ) 2 + (z - 24) 2 = d\ 

• KNN: For this algorithm, we considered the 
value of k = 2 and 3. Where k represents the 
number of nearest neighbors. In addition, 
Euclidean distance was used as a measure of 
similarity; 

• Artificial neural networks (MLP and RBF): 

Table. 3: Parameters Adopted for the Neural Networks 


Parameters 

Value 

Number of neurons in the 
input layer 

20 

Number of neurons in the 
output layer 

3 

Number of neurons in the 
hidden layer 

50,100.150 
200.250 r 300 

Number of hidden layers 

2 MLPe1 RBF 

Activation 

Function 

Sigmoide 

Training algorithm 

Back propagation 

Learning rate 


Number of iteration 

1 5000 


Probabilistic Model (Pure version of Bayes theorem): In 


using the histogram and kernel based versions. The kernel 
method was modeled by a distribution by a normal 
distribution, according to equation (6). 

4.3 RESULTS OF SCENARIO I 

For this scenario, it was considered 150 RP's distributed 
in the two floors. For each PR, 300 RSS values were 
collected, according to the following rule: Notebook with 
heights of 0.2m; 1.0m and 1.80m. In each height, 25 
measurements were performed with a notebook with the 
Notebook facing toward the north, east, south, and west 
directions. The Vistumbler software [25] was used to 
acquire RSS. 

After reading the RSS in all RPs, an experiment was 
conducted to verify the variation of the signal in the 
environment. At a fixed point (15.2.0.7 and 5.72m), for a 
period of 300s the signal strength Signal strength was 
observed with respect to the four APs. The results are 
shown in the figure 6. Note that even from a fixed 
position, the fluctuations in the signal level reaches 10 
dB, which is the expected variability for indoor 
environments. 

With respect to the performance of the algorithms, the 
method of lateration when compared with others was not 
efficient. This is an expected result, since methods based 
on traditional outdoor location are not suitable for indoor 
environments, which have a much greater variability [11]. 
The system implemented with this method showed a 
mean error e m of 2.23m from the actual locationwith 
precision of 68 % and 80.1% relative toabsolute 
maximum errorse a of 1.5 and 3.5m respectively. 


API: (X = 8.5, Y = 1.72, 2 s 3.69) 


| no 

T3 

c 
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* I.' 
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J l|V|I v UI W 

I I 1 | '» 

f IJ 

v| 


i i 
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(a) 


AP2: (X = 12.8, Y = 1.5. Z = 8.0) 



AP3:(X = 7.45. Y = 1.6, Z = 6.93) 



AP4: (X = 13.1. Y = 1.7, Z = 5.07) 



<d) 


Fig. 6: Variations of RSS values from four APs in a fixed point 
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Table.4: Statistics of positioning estimates 


ALGORITHM 

«. <*> 

^ 133m 

(%) 

^ 3.5 rn 

^ 5.5m 

<*'(%) 

« 7.5m 

Lateration 

68,2 

80,1 

87,8 

95.1 

kNN s 2 

75,8 

88 

91,8 

98.9 

kNN s 3 

76.1 

88.7 

95.3 

99.1 

RNA-MLP 

76,2 

89 

95,8 

99.3 

RNA•RBF 

76,65 

89.6 

96,5 

99,4 

HISTOGRAM 

76,9 

90.8 

98 

100 

KERNEL 

77,2 

91,2 

98,6 

100 

PROPOSED 

METHOD 

77,9 

93.1 

98,9 

100 


Table. 5: Statistics of positioning estimates (average, 
coefficient of variation, Ql, Q2, and Q3 quartiles in 
meters) 


ALGORITHM 


C. 

Ql 

Q2 

Q3 

Lateration 

2.235 

59.70% 

1.36 

1.73 

2.57 

kNN = 2 

1.919 

15.22% 

1.33 

1.66 

1.93 

kNN = 3 

1.916 

41.10% 

1.33 

1.66 

131 

RNA-MLP 

1.905 

43.15% 

1.33 

1.00 

1.90 

RNA-RBF 

1.883 

41.88% 

1.33 

1.65 

1.88 

HISTOGRAM 

1.813 

37.97% 

1 32 

1 65 

~TS0 

KERNEL 

1.832 

30.88% 

1.32 

1.65 

1.81 

PROPOSED 

METHOD 

t.803 

35.06% 

1.32 

1.64 

1.81 


0.95 


| 0.9 

15 


0.85 


I 

3 0.8 

I 

3 

o 

0.75 


0.7 




1 R l \NC,l 1 \TION 

-m- 

KNN *2 


KNN * 3 


ANN - MLP 


ANN - RBF 

- # 

BAYES - HISTOGRAM 

- • 

BAYES - KERNEL 

- » - 

PROPOSED METHOD 


1.5 3,5 5,5 7,5 9,5 11,5 

Distance error (Meters) 

Fig. 7: Error in meters - cumulative 


13,5 


estimation of the RSS attenuation applied to the lateration 
method, this simulation includes the estimation of the 
vectors RSS in each PR. It was observed that in the 
scenario in question, the lateration method presented 
performance superior to that implemented in scenario I. 
This better performance can be explained by the smaller 
variation in the signal compared to scenario I, because 
although the propagation model used in this work is 
considered more complete with respect to the others, in 
order to include the losses from the various obstacles, can 
not accurately represent the signal fluctuations in the 
various types of indoors environments. Even with better 
performance in this scenario, the lateration method was 
less accurate than the others, presenting e m of 2.0m, 
against the most similar algorithm in terms of 
performance (kNN with k = 2) that it obtained e m of 
1.92m. In relation to the other algorithms, there were no 
significant differences in performance (Algorithms based 
on the fingerprint method, are less influenced by the 
variation of RSS in the environment). 
The proposed method presented^ equal to 1.786m. The 
results are summarized in the tables 6, 7 and in the figure 
8 . 


Table.6: Statistics of positioning estimates (average, 
coefficient of variation, Ql, Q2, and Q3 quartiles in 
meters) 


ALGORITHM 

1 «.(%r 

^ 1.5m 

<»(%) 

^ 3.5m 

C, (%) 

^ 5.5m 

<« (%) 

^ 7.5m 

Lateration 

72,9 

87,4 

93,70 

97,8 

kNN -2 

76 

88,1 

95 

99 

kNN = 3 

76,05 

88,5 

94,98 

99 

RNA-MLP 

76,5 

88,7 

95,55 

99,7 

RNA-RBF 

76,66 

89,5 

96,3 

99,8 

HISTOGRAM 

77 

91 

97,9 

100 

KERNEL 

77,1 

91,3 

98.5 

100 

PROPOSED 

METHOD 

79,1 

93 

98,8 

100 


4.4 RESULTS OF SCENARIO II 

In this scenario, equation (13) was used to simulate the 
signal behavior in the environment. In addition to the 
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Table.7: Statistics of positioning estimates (average, 
coefficient of variation, Ql, Q2, and Q3 quartiles in 
meters) 


ALGORITHM 


c v 

Ql 

Q2 

Q3 

Lateration 


5U?0%~ 

1 34 

1.69 

2.00 

kNN = 2 

1.921 

44.89% 

1.33 

1.66 

1.92 

kNN = 3 

1.908 

44.92% 

1.33 

1.66 

1.91 

RNA - MLP 

1.897 

13.59% 

1.33 

1.66 

1.92 

RNA-RBF 

1.880 

12.38% 

1.33 

1.65 

1.88 

HISTOGRAM 

LSI! 

38.13% 

1 32 

1.65 

1.85 

KERNEL 

1.791 

35.07% 

1.31 

1.63 

2.00 

PROPOSED 

METHOD 

1.7*6 

34.54% 

1.31 

1.63 

2.00 



1.5 3,5 5,5 7,5 9,5 11,5 13,5 

Distance error (Meters) 

Fig..8: Error in meters - cumulative 


V. CONCLUSION 

In this paper, we propose a three-dimensional (3D) 
solution to this problem through a hybrid algorithm. In 
addition, a comparison was made with the main 
algorithms discussed in previous research papers, 
including with the lateration method, one of the most used 
in the outdoor context. With the exception of this method, 
all algorithms used a fingerprint recognition technique. 

A lower variation in signal levels in Scenario II had a 
direct influence on the performance of the lateration 
method. In this scenario, the lateration obtained better 
performance than in scenario I (around 23cm in relation 
to e m ). Even with better performance in this scenario, the 
lateration method was less accurate than the others, which 
also demonstrated better performance compared to 
scenario I, but not significantly. As discussed, indoors 
environments have high variability. IPS developed only 
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from simulations based on propagation models may not 
capture all the information regarding the signal variations 
in this type of environment and therefore, overestimate 
the results, as in the case of the lateration method 
implemented in this work. 

The proposed algorithm exhibited a better performance 
than that of other algorithms, resulting em = 1.80 and 
1789 for scenarios I and II respectively. Because the 
proposed algorithm is based on Bayesian inference, there 
is good integration with the fingerprint technique, this 
because it relates the current and past information 
efficiently, which results in good results in relation to the 
classification. However, one of the limitations of the 
techniques in question is related to the quality of the 
information, that is, a robust database is necessary for an 
IPS based in this algorithm to present good results. 

As future work, we intend to develop the following 
solutions in order to maximize the precision of the 
location estimation: 

• IPS based on angulation; 

• Hybrid IPS based on triangulation and 
fingerprint; 

• IPS based on multiple discriminant analysis. 
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